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Abstract 

There are two mechanisms for the generation of an asymmetry between the strange and anti- 
strange quark distributions in the nucleon: nonperturbative contributions originating from nucleons 
fluctuating into virtual baryon-meson pairs such as AK and T,K, and perturbative contributions 
arising from gluons splitting into strange and anti-strange quark pairs. While the nonperturbative 
contributions are dominant in the large- x region, the perturbative contributions are more significant 
in the small- x region. We calculate this asymmetry taking into account both nonperturbative and 
perturbative contributions, thus giving a more accurate evaluation of this asymmetry over the 
whole domain of x. We find that the perturbative contributions are generally a few times larger in 
magnitude than the nonperturbative contributions, which suggests that the best region to detect 
this asymmetry experimentally is in the region 0.02 < x < 0.03. We find that the asymmetry 
may have more than one node, which is an effect that should be taken into account, e.g. for 
parameterizations of the strange and anti-strange quark distributions used in global analysis of 
parton distributions. 
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I. INTRODUCTION 



Partons in the sea of the nucleoli play an important role in understanding many properties 
of the nucleon and in explaining many experimental results in hadronic physics. For example, 
it is now well established [IH1] that gluons carry about 1/2 of the nucleon momentum, while 
the valence and sea quarks carry the other half. It is also well known that about 1/3 
of the nucleon spin can be attributed to the spin of the valence quarks, while the rest 
should be attributed to the spin and orbital angular momenta of quarks and gluons [5, 6J. 
Also, subprocesses involving strange and charm quarks of the nucleon contribute 20% to W 
production at the LHC [7J. A precise understanding of the strange and charm content of 
the nucleon is important to the search for physics beyond the Standard Model. 

The asymmetry between the strange and anti-strange quarks of the nucleon (s-s asym- 
metry) [8] could affect the extraction of the Weinberg angle from neutrino-nucleon deep 
inelastic scattering [9] . Most theoretical estimations for this asymmetry using quark models 
of the nucleon suggest (see e.g. [TO]) that this effect alone would not be large enough to 
explain the discrepancy between the NuTeV measurement of sin 2 6\y [IT] and the current 
world-average value p2]. Global analysis of hard-scattering data generally provides weak 
constraints on the s-s asymmetry of the nucleon [13-23J, mainly due to the fact that there 
is limited data sensitive to the strange content of the nucleon. Some groups have assumed 
symmetric strange distributions, i.e. s = s, in their latest global fits for the parton dis- 
tribution functions (PDFs) of the nucleon [TBTflE] , while other groups were able to find a 
strangeness asymmetry that has the magnitude large enough to remove the NuTeV anomaly 
(see e.g. [2TH25] ) . In determining the strangeness asymmetry, most groups adopted func- 
tions with one node as the nonperturbative input for the asymmetry at the initial scale 
[T3] [HI [TH [201 121], although there is no physical principles preventing the asymmetry hav- 
ing multiple nodes. On this aspect, the neural network technique employed by the NNPDF 
group has the advantage of allowing for multiple nodes [2~TH2"3"] . 

We can classify sea partons according to the dynamics responsible for their creation: 
perturbative contributions originating from gluons splitting into quark-antiquark pairs, and 
nonperturbative contributions, which come from nucleons fluctuating into baryon-meson 
pairs with the partons in the baryon-meson pairs manifesting as the sea partons of the nu- 
cleon. The perturbative contributions can be investigated using QCD evolution equations for 
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the parton distribution functions (PDFs) of the nucleon; the nonperturbative contributions 
can be estimated using nonperturbative models for the nucleon. At leading order and next 
to leading order of the strong coupling constant a s , the perturbative contributions usually 
respect symmetries of quark models of the nucleon, such as flavor symmetry, quark- ant iquark 
symmetry and charge symmetry, whereas the nonperturbative contributions generally break 
these symmetries [26] . 

The breaking of the symmetry between the strange and anti-strange quark distributions 
of the nucleon was predicted over two decades ago by Signal and Thomas [E] using the 
meson cloud model (MCM) of the nucleon [27] . The asymmetry they investigated has the 
nature of a nonperturbative contribution. Many theoretical papers have since investigated 
these nonperturbative contributions 0, HOI EHHS] • These calculations generally predict the 
distribution x(s — s) peaks in the region of 0.1 < x < 0.2. 

Although it has been generally assumed that perturbatively generated sea distributions 
respect the symmetries of quark models of the nucleon, it was shown recently [I2J that 
the quark-antiquark symmetry of the sea is violated when the perturbative QCD evolution 
is calculated at next-to-next-to-leading order (NNLO). This occurs because the splitting 
functions for quarks and antiquarks are different from each other at that order [43] . To 
illustrate this evolution effect, a null symmetry at the initial scale, i.e. (s — s)(x,Qq) = 0, 
was assumed in [42J. The calculated distribution x(s — s) shows a maximum at x < 0.1. 

In this work we calculate the strange and anti-strange asymmetry in the whole x region, 
including both perturbative and nonperturbative contributions to the asymmetry. We briefly 
summarize the formalism for calculating the perturbative and nonperturbative contributions 
to the asymmetry in Sec. 2. The numerical results and discussion, including possible 
implications for the NuTeV anomaly, are presented in Sec. 3. The last section is reserved 
for a summary. 

II. FORMALISM 

A. Nonperturbative Contributions 

In the meson cloud model the nucleon can be viewed as a baryon 'core' surrounded by 
a mesonic cloud. The wave function of the nucleon can be expanded in terms of a series of 
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baryon plus meson Fock states 

N) phys = Z\N) haTe + J2J2 ! dyd 2 k ± ^ B x M (y,k 2 ± )\BM), (1) 

BM XX' J 

where Z is the wave function renormalization constant, (p B x ^(y, k\) is the wave function 
of the Fock state containing a virtual baryon (B) with longitudinal momentum fraction y, 
transverse momentum k^, and helicity A, and a virtual meson (M) with momentum fraction 
1 — y, transverse momentum — kj_, and helicity A'. 

The lifetime of the virtual baryon-meson components is generally much longer than the 
interaction time in the deep inelastic process, thus the quarks and antiquarks in the baryon 
and meson contribute to the parton distributions of the nucleon. These nonperturbative 
contributions can be expressed as a convolution of fluctuation functions with the valence 
parton distributions in the baryon B or meson M. The contributions to the asymmetry 
distribution at the scale Q, (s — s)(x,Q 2 ), can be written as, 

^dy 

y 



(s - s)(x, Q 2 



x 



fBM/N(y)s B (-,Q 
y 



-fMB/N(y)s M (-,Q 2 
y 



(2) 



where B = A(S), M = K(K*), sb and % represent the distributions of strange quark and 
anti-strange quark in the baryon and meson, respectively. The f B M/N(y) (/ms/at(?/)) is the 
fluctuation function which describes the probability of a nucleon fluctuating into a baryon 
(meson) with longitudinal momentum fraction y. We have f bm /n(]J) = Im b /n ( 1 — y) due 
to the conservation of momentum and charge, 

roc 

fBM/ N (y) = J2 / dk i<f>BM(y, k 2 ± W B M(y, k 2 ± ). (3) 

AA' J ° 

The baryon-meson wave function, k 2 ^), and thereby the fluctuation functions 

can be calculated using effective meson-baryon-nucleon interaction Lagrangians 



£nbp = igNBpNj 5 PB, (4) 
Cnbv = g N BvN lf y»B + f NBV Na, u B(d^ - d v V»), (5) 

where gNBP and Qnbv are the coupling constants. iV and B are spin 1/2 fields, P and V 
are pseudoscalar and vector fields, respectively. The anti-symmetric tensor a^ u is defined as 



°V = *[7/xj7f]/2- For the latest calculations for the asymmetry (s — s) using the MCM and 
more calculation details we refer the reader to [TU] . 

B. Perturbative Contributions 

The perturbative contributions to the strange and anti-strange quarks in the nucleon 
arising from gluons splitting into quark-antiquark pairs can be taken into account using 
perturbative evolution in QCD. Up to next-to-leading order (NLO) in a s the probability of a 
splitting q — > q' is identical to that of q — > q'. A difference between the probabilities for these 
splittings arises at NNLO, which consequently leads to a strange-anti-strange asymmetry 



The splitting function P a ^ describing the splitting b — > a is expanded in terms of a 
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,4:71 
71=1 V 

The splitting functions have been calculated up to order n = 3 |33j |33] . Under the assump- 
tion of charge conjugation invariance and flavor symmetry, the splitting functions P a b can 
be written as [331 133], 

Pqiqk = Pq~iq~k = ^ikPqq + Pqq j (7) 
-^"Wfc = Pq~i<lk = ^ikPqq + Pqq ■ (8) 

where the functions P qq and describe splittings in which the flavor of the quark changes. 
The relation P qq = P qq holds up to NLO in a s . 

Only the evolution of flavor nonsinglet (NS) combinations of parton distributions needs 
to be considered in order to evaluate (s — s) since it is a flavor NS quantity. The three NS 
distributions considered in [42j are 

n f 

/ v = £(/«-/*)> ( 9 ) 

i=i 

1 n/ 

fqi = hi ^ fqi ~ ~ ^ (/<?;) i f< 

" i=i 

where rif is the number of flavors. 

The evolution equations for these NS distributions are 
d f j {x,Q 2 ) f 1 dz 



,„~ ■ (10) 



rflnQ 2 



/ d -P>\ X -,a s {Q*)\ f(z,Q 2 ) (11) 
j x z iz J 



with j = ±, V. The evolution kernels appearing in Eq. (11) are given by 

P ± 



p v ± p v 

m m' 



qq 



qqj 



p: 



Air 



p(2)S 



(12) 



where Pns is the three-loop non-singlet splitting function. 

It is convenient to solve the evolution equations in Mellin space where the moments of 
the distributions f^{Q 2 ) and the splitting functions P^(x), 



dxx N - 1 f j {x,Q 2 



pi 

1 x — 



dxx^PHx) 



(13) 
(14) 



are used. 

The Mellin moments of the distribution (s — s) were found 



to evolve according to, 



[s - s) N (Q 2 ) = U-(Q,Q )(s-s) N (Ql) 



P. 



(2)5 



ns,N 



2f3 N f 



47T 



<*s(Ql) 

An 



x (u v + d v ) N (Q 2 ), (15) 

where the notations f s — s and f§ = s have been used, and u v and d v represent valence 
w-quark and <i-quark distributions respectively. 



The first term in Eq. (15) represents nonperturbative contributions to the asymmetry 



while the second term represents the perturbative contributions. The evolution operator 

" 2 dq 2 



U (Q, Qo) = exp 



-PnsMq 2 )} 



(16) 



evolves the nonperturbative asymmetry at Qq to higher values of Q 2 . 

In Ref. jl2] the assumption (s — s)jv(Qq) = was adopted and thus only the perturbative 
contributions to the asymmetry were studied. The predictions for the asymmetry in the 
x-space, (s — s)(x, Q 2 ) were obtained via a numerical Mellin inversion. 



III. TOTAL STRANGENESS ASYMMETRY 



In order to evaluate the asymmetry in the whole region of x we need to take into account 
both perturbative and nonperturbative contributions. Rather than using a numerical Mellin 
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inversion, we apply the inverse Mellin transformation to Eq. (15) to obtain an expression in 
x space for the total asymmetry, 



(s - sf°\x, Q 2 ) = (s - s) NP (x, Q 2 ) + (s- s) p (x, Q 2 ), (17) 

with 

X f^*(;) (u V + d V ) (z,Q 2 ). 

(18) 

( 2) S 

Both the exact expression and a parameterization form for Pns (x) are given in [33]. The 
parameterization form, which deviates from the exact expression by less than one part in 
thousand, is sufficiently accurate for our calculations, 

pV> s (x) ^ n f {[Loi-lWAx- 1 - 7.208x) + 151.49 
+44.51x - 43. 12x 2 + 4.82x 3 ] [1 - x] 
+L L 1 [-173.1 + 46.18L ] + 178.04L 
+6.892L 2 + 40/27[L^ - 2L^\ } , (19) 

with Lq = Inx and L\ — ln(l — x). 

We use the result obtained in the MCM (Eq. ([2])) to estimate the nonperturbative con- 
tributions, (s — s) NP (x, Q 2 ). 

IV. NUMERICAL RESULTS AND DISCUSSION 

In the numerical calculations for the nonperturbative contributions we considered Fock 
states AK, AK*, Y,K, and UK*, and used the parameters given in [10]. The strange and 
anti-strange quark distributions in the baryons and mesons were calculated using the MIT 
bag model developed by the Adelaide group [45TT4T] and Massey group [IB] . The distributions 
were evaluated at an initial scale of fi^ = 0.23 GeV 2 and evolved to Q 2 = 16 GeV 2 using the 
program given in [49]. 

The results, in comparison with the asymmetry obtained with the NNPDF2.3 PDF set 
[23] . are shown in Fig. [TJ One can see that the non-perturbation calculations for the asym- 
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FIG. 1: Nonperturbative contributions to the asymmetry x(s — s) at Q 2 = 16 GeV 2 calculated with 
the meson cloud model, in comparison with the total asymmetry obtained with the NNPDF2.3 
NNLO PDF set. The solid and dashed curves are the MCM results considering iT-meson contri- 
butions and K- and i^*-mesons contributions, respectively. The shaded area represents one-sigma 
uncertainty band for the NNPDF2.3 set. 



metry are generally much smaller than the total asymmetry obtained with the NNPDF2.3 
PDF set. 



The perturbative contributions depend on the initial scale chosen (see Eq. (18)). Two val- 
ues for the initial scale, Qo = 0.51 GeV and Qq = 1.1, were used in our numerical evaluations 
as in [12] . It is worth noting that the calculations depend on the up and down valence dis- 
tributions at high Q 2 where hard scattering data available, rather than the distributions at 
the initial scale Qq. The valance distributions of the nucleon are reasonably well determined 
via global fits. Although the d v and u v obtained with different PDF sets may not agree 
within the quoted uncertainties, the differences are generally smaller than 10%. Thus we 
can expect similar level of uncertainties due to the choice of PDF sets in the calculations for 



the perturbative contributions to the strangeness asymmetry using Eq. (18). We employed 
the ABM11 NNLO PDF set for the d v and u v at Q 2 = 16 GeV 2 in our evaluations. 

The perturbative contributions to the asymmetry x(s — s), in comparison with the asym- 
metry obtained with the NNPDF2.3 PDF set [23], are shown in Fig. [2j The results are in 
agreement with those presented in [12]. Note that we have used an inverse Mellin transfor- 



mation of Eq. (15) to obtain the expression for perturbative contributions in the x space 
(Eq. (18)) while the Mellin inversion was done numerically in [12]. The perturbative con- 



tributions in the small x-region (x < 0.05) tend to agree with results obtain with the 
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FIG. 2: Perturbative contributions to the asymmetry x(s — s) at Q 2 = 16 GeV 2 , in comparison 
with the total asymmetry obtained with the NNPDF2.3 NNLO PDF set. The solid and dashed 
curves are the results calculated with Qq = 0.51 GeV and Qq = 1.1 GeV, respectively. The bands 
in the perturbative calculations represent the one-sigma range due to the uncertainties associated 
with the d v and u v of the ABM11 PDF set. The shaded area represents one-sigma uncertainty 
band for the NNPDF2.3 set. 

NNPDF2.3 NNLO PDF set, which suggests that the perturbative contributions are the 
dominant mechanism for the generation of this asymmetry at small x. 

From Figs. [T] and [2] one can see that the nonperturbative contributions dominate in the 
region x > 0.1 while the perturbative contributions dominate in the region x < 0.1. The 
nonperturbative contributions are smaller in magnitude than the perturbative contributions, 
which indicates that the asymmetry may be most easily measured in the region of x around 
0.03. 

The calculations for the total asymmetry including both perturbative and nonpertur- 
bative contributions are shown in Figs. [3] and |4j Only the i^-meson contributions were 
considered in the nonperturbative calculations for Fig. |3j while both K- and iP-mesons 
contributions were included in the nonperturbative calculations for Fig. |4| 

Figs. [3] and [4] suggest that the asymmetry may have more than one node, in contrast to 
the assumption of only one node for the asymmetry commonly used in most global analyses 
of the PDFs [131 IH EE ED]- In the case of one node, the calculated asymmetry presented 
in this work can be described using the parameterization suggested in 



s_ ( x , Q) = A x Al (1 - x) A2 tan 



cx 



X 



dx+ex 



(20) 



where A , A\, A 2 , a, b c, d, and e are parameters. It is very difficult to describe our calcu- 
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FIG. 3: The total asymmetry x(s — s) at Q 2 = 16 GeV 2 including both the perturbative contri- 
butions and nonperturbative contributions which are calculated with only K mesonic cloud. The 
solid and dashed curves are the results obtained with Qq = 0.51 GeV and Qq = 1.1 GeV in the 
perturbative calculations. The bands represent the one-sigma range due to the uncertainties asso- 
ciated with the d v and u v at Q 2 = 16 GeV 2 . The shaded area represents one-sigma uncertainty 
band for the NNPDF2.3 set. 
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FIG. 4: Same as Fig. [3] but both K and K* mesonic clouds are considered in the nonperturbative 
contributions. 



lated asymmetry using other forms simpler than Eq. (20) and involving fewer parameters 



It is interesting to note that when only i^-mesons are considered in the nonperturbative 
contributions the calculated asymmetry exhibits a shape similar to that determined with 
the NNPDF2.3 NNLO PDF set (see Fig. [3j) , though the magnitudes at the region of x > 0.2 
are rather different. 

The current experimental constraints for the strangeness of the nucleon come mostly from 
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TABLE I: (x(s - a)) at Q 2 = 16 GeV 2 



Qo(GeV) 


P 


P+NP without K* 


P+NP including K* 


0.51 


-5.39 x 10~ 4 


-4.09 x 10~ 4 


-7.09 x 10~ 4 


1.10 


-1.27 x 10~ 4 


3.33 x 10~ 6 


-2.96 x 10~ 4 



the neutrino dimuon production data, and are not strong enough to allow a meaningful fit 



for the strangeness distributions using the many numbers of parameters of Eq. (20) jT 
This highlights the need for more accurate future experiments to constrain the strange and 
anti-strange distributions of the nucleon. The Large Hadron Collider has the great potential 
to provide new information for the strange content of the nucleon, as shown by a recent 
study by the ATLAS Collaboration [50] . 

To investigate the possible effects on the NuTeV measurement of the Weinberg angle 
we calculated the second moment of the asymmetry (x(s — s)} = dxx(s — s). The 
results are given in Table |T] where P and NP stand for the perturbative and nonperturbative 
contributions. The results depend on the model parameters - the initial scale chosen in the 
perturbative calculations and whether the K* mesonic cloud is included in the calculations 
for the nonperturbative contributions. A value of (x(s — s)) ~ 0.004 would be required if the 
anomalous NuTeV result was caused only by the strange quark asymmetry. Our calculations 
suggest that the second moment (x(s — s)) is negative, but is about an order of magnitude 
too small to have any significant effect on the NuTeV result for the Weinberg angle. Thus 
other mechanisms, as discussed in [HI I5TH53] , could be responsible for the NuTeV anomaly. 



V. SUMMARY 



In this paper, we analyzed the strange and anti-strange quark asymmetry in the nucleon 
sea. We have taken into account both the nonperturbative contributions as well as the con- 
tributions from perturbative QCD evolution. The nonperturbative contributions were cal- 
culated using the meson cloud model in which the baryon-meson Fock state, A" — > AK(K*), 
and A" — > HK(K*) were included. The perturbative contributions arises at the NNLO in 
perturbative QCD because the splitting function for the splitting q — > q' differs from that 
for the splitting q — > q' starting from that order. 

The nonperturbative contributions are dominant in the region x ^ 0.1 while the pertur- 
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bative contributions are dominant in the smaller x region. It was found that the perturbative 
contributions are generally a few times larger in magnitude than the nonperturbative con- 
tributions, which suggests that the best region to detect this asymmetry experimentally is 
in the region of x around 0.03. We found that the asymmetry may have more than one 
node which could shed light on the parameterization of the strange and anti-strange quark 
distributions in the global analysis of the parton distribution functions. More experiments 
directly detecting the strangeness of the nucleon are highly desired. We found that the 
asymmetry is not large enough to have any significant effects on the NuTeV measurement 
for the Weinberg angle. 
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